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SUMMARY 


The near-24-hour Syncomll satellite (with an almost circular orbit) has been under continuous 
observation by range and range-rate radar and minitrack stations for 7 months since mid- August 
1963, when the orbit was relocated, placing its mean longitude at about 55 degrees west of Green- 
wich. During the first 4 months of this period, the satellite was allowed to drift free in the grav- 
ity fields of the earth, sun, and moon. In this first free-drift period, the satellite experienced a 
mean daily drift acceleration of its ascending node (with respect to Greenwich) of 

- (1.27 ± .02) x 10~ 3 degrees/day 2 . (1) 


The average growth of the semimajor axis for this period was 

(.0993 ± .0042) km/day . (2) 


These values, checked by a simulated particle trajectory run on the Goddard ITEM program, con- 
firm a significant longitude-dependent earth-gravity potential. The existence of a "triaxial earth" 
has been a speculation of geodesists since the early years of this century. 


During the last 3 months of this 7-month drift period, starting at the end of November 1963, 
Syncom n was relocated atiibout^&O degrees west of Greenwich. In this period, the mean daily 
drift acceleration of the a^cendiiig ripde was 


v 


- (1.32 ± .02) x 10 3 degrees/day 2 


(3) 


The average growth of the semimajor axis for this period was 


(.0994 ± .0080) km/day . 


(4) 


Combining the results of (1) through (4) above for the two separate drift periods, it is estimated 
(on the basis of a triaxial geoid only) that the absolute magnitude of the longitude dependent- gravity 
coefficient J 22 is 

J 22 = - (1.7 ± .05) x 10~ 6 . 


This value corresponds to a difference in major and minor earth- equatorial radii of 


(a c - b c ) = 213 ± 6 feet . 


The best present estimate of the position of the earth's major equatorial axis is 

k 22 = 19 ± 6 degrees west of Greenwich . 
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In view of previous estimates of the higher order tesseral harmonics of the earth's field, the 
true value of J 22 , separated from the small influence of gravity anomalies of third and higher 
order on the reduction for a triaxial earth only at "synchronous" altitudes, will probably be some- 
what higher than the -1.7 x 10' 6 reported herein. The true value of J 22 , however, is not likely to 
be greater than -2.2 x 10~ 6 or smaller than -1.6 x 10" 6 . The true location of the earth equator's 
major axis is not expected to differ significantly from that reported herein, when all higher tes- 
seral harmonics are accounted for. (See Appendix B.) 

The reported value of J 22 = -1.7 x 10' 6 implies that a maximum tangential velocity correc- 
tion of 

av t = 5.36 ft./sec./year 

is required to keep a satellite with a 24-hour circular equatorial orbit continuously "on station" 
at a longitude midway between the longitudes of the equatorial major and minor axes of the earth. 
The original "conservative” Syncom I design requirement of av t = 17 ft./sec./year correction 
capability was based on the longitude -dependent earth field of Izsak (January 1961), which is now 
outdated. (See Appendix B and Reference 4.) 


ii 



CONTENTS 


Section Page 

INTRODUCTION: SUMMARY OF PREVIOUS INVESTIGATIONS AND 

DISCUSSION OF RESULTS 1 

1. THE BASIC THEORY OF THIS REDUCTION 2 

2. EVALUATION OF THE PERTURBING FORCE 5 

3. COMPLETION OF THE DERIVATION OF THE DRIFT EQUATIONS 8 

4. GENERAL CONSIDERATIONS OF THE SOLUTIONS OF THE DRIFT EQUATIONS . . 10 

5. APPROXIMATIONS TO THE EXACT DRIFT SOLUTIONS FOR PERIODS VERY 

CLOSE TO SYNCHRONOUS 13 

6. DETERMINATION OF EARTH- EQUATORIAL ELLIPTICITY FROM TWO 

OBSERVATIONS OF DRIFT ACCELERATION AT A GIVEN LONGITUDE 
SEPARATION 16 

7. REDUCTION OF 27 SYNCOM H ORBITS TO DETERMINE THE EARTH'S 

EQUATORIAL ELLIPTICITY 18 

REFERENCES 25 

LIST OF SYMBOLS 27 

Appendix A— Reduction of Simulated Particle Trajectories for Earth- Equatorial 

EUipticity 29 

Appendix B— The Earth Gravity Potential and Force Field Used in This Report: 

Comparison with Previous Investigations 35 

Appendix C— Expressions for The Inclination Factor 41 

Appendix D— Derivation of The Exact EUiptic Integral of Drift Motion for A 24- Hour 
Near-Circular Satellite Orbit: Comparison of The Exact Solution With 
The Approximate Solutions For Periods Very Close To Synchronous 45 

Appendix E— Basic Orbit Data Used in This Report 49 

Appendix F— Calculation Of The Radiation Pressure On Syncom H 51 

iii 



INTRODUCTION 


SUMMARY OF PREVIOUS INVESTIGATIONS AND DISCUSSION OF RESULTS 

The question of the existence and extent of the longitude- dependence of the earth's external 
gravity field has concerned geodesists since the early years of this century. (See Reference 1 for 
example.) The existence of a longitude- dependent field implies the existence of inhomogeneities 
and states of stress within the earth which are of fundamental importance to all dynamical theories 
of the earth's interior. 

In Appendix B, I have summarized 23 reductions of gravimetric, astro- geodetic, and satellite 
gravity data reporting longitude-dependent terms in the earth’s external gravity field represented 
as a series of spherical harmonics. Ten of these reductions are based on worldwide surface- 
gravity measurements only. Although surface measurements have the advantage of providing an 
excellent sampling of the field in latitude and longitude, they are seriously affected by even small 
uncertainties in station position with respect to the mean earth geoid, as these are of the same 
order- of- magnitude as the reported geoid deviations caused by longitude- dependent gravity. 

Of the 12 satellite reduced tesseral gravity fields, those of Kozai (1962), Izsak (July 1963), 
Kaula (September 1963), and Guier (1963) show good general agreement in tne eastern hemisphere 
when the constants are used as a set. All of these observers are aware, nevertheless, of the high 
degree of uncertainty in the reported values of the individual coefficients themselves, this un- 
certainty being due mainly to unresolved station- datum errors and to the limited sampling of the 
field from observations on a small number of medium- altitude, medium- inclination satellites. 

As late as July 1963, Izsak stated (Reference 2): "It might be some time before one can ar- 
rive at definite conclusions regarding the longitude dependence of the earth's gravity field.” The 
presence of Syncom n, high over the earth with a nearly stationary, narrow figure- 8 ground track 
centered close to the equator, dramatically alters this gloomy outlook. The 24-hour satellite is 
high enough to be unaffected by the earth's atmosphere, yet is close enough to the earth to be pro- 
tected from the solar wind by the "magnetosphere" of the earth, and to remain essentially unaf- 
fected by sun or moon gravity. 

In theory, as the ground track is nearly stationary, any small earth-gravity anomaly in longi- 
tude will in time, cause significant drift of the ground-track configuration. In theory too, only 
observations of the longitude of the satellite from a single ground station are necessary to reveal 
this effect of the "tesseral" gravity field over an extended period of time. Long-term observa- 
tions of the longitude drift of one or more 24-hour satellites should reveal the exact nature of the 
tesseral- gravity field to at least the third order without essential difficulty. The great height of 
the 24- hour satellite tends to cancel out the individual contributions to the longitudinal drift of 
anomalies higher than about the fourth order. It is fortunate that the initial slow westward drift 
of the ground track of Syncom n (August 1963 to March 1964) appears to have occurred relatively 
close to a point midway between the triaxial earth's major and minor equatorial axes where the 
perturbation of the second tesseral-gravity anomaly, for which the reduction was made, is great- 
est. A weighted average of the longitude-perturbation fields of Kozai (1962), Zhongolovitch (1957), 
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Kaula (September 1963), and Izsak (July 1963), at the altitude and longitude of Syncom II during 
this period, shows that the second tesseral should be contributing about 80 to 85 percent of the 
perturbing force. If further observations of this and other 24- hour satellites confirm this esti- 
mate, then the magnitude of the dominant J 22 potential term, separated from higher order effects, 
will increase by 25 percent at most from the value reported here, which is based on the assump- 
tion that only the triaxial earth- gravity field is being measured by the drift observations. The 
reported location of the major equatorial axis of the triaxial earth ellipsoid is not expected to 
change significantly with this later refinement. 

In summary, long-term observations on the drift of Syncom n have already established to a 
high degree of certainty that: 

1. The earth must be considered to be a "triaxial ellipsoid" (for example, having a sea-level 
surface of this form) for the purposes of 24-hour satellite design. (For broader geodetic 
purposes; a significant longitude dependent gravity field exists, defining: 

2.2 x 1(T 6 < |j 22 1 < 1.6 x 1(T 6 ; -25° < ^ 22 < -13°.) 

2. The difference between the major and minor equatorial radii of that ellipsoid is not less 
than 200 feet nor greater than 225 feet. 

3. The location of the major equatorial axis of the "triaxial geoid" is between 13 degrees and 
25 degrees west of Greenwich. 

It may be added that the study of simulated 24-hour satellite drift in a triaxial earth field, influ- 
enced also by sun and moon gravity and by sun- radiation pressure perturbations, shows that the 
theory of longitude drift presented in this report is substantially unaffected by all perturbations 
except that due to the earth's elliptical equator and possible higher order longitude-dependent 
earth-gravity anomalies. 

1. BASIC THEORY OF THIS REDUCTION 

(Determination of the Longitude Drift and Orbit Expansion for a 24-Hour Satellite With a 
Near-Circular Orbit Affected by a Small But Persistent Tangential Pejtdrbing Force.) 

The dominant perturbations of a 24- hour equatorial satellite in a higher order earth- gravity 
field have been derived many times in the literature (References 3, 4, 5, and 6). In these refer- 
ences, the perturbations were found by directly linearizing the equations of motion themselves and 
displaying the perturbed motion in appropriate geographic coordinates; no attempt was made to 
treat the drift of the inclined 24-hour satellite. 

In this report, I will depart from this rather involved and difficult- to- visualize procedure of 
linearization-of-the-equations-of-motion. Instead, I will show how simple it is to derive the 
dominant drift and orbit- expansion equations for the 24-hour satellite by dealing with what can be 
called "the perturbation of the 2-body energy" of the geographically stationary satellite by the 
small but persistent longitude -dependent earth-gravity force. This paper will not discuss in 
detail the limits of validity of the expressions derived. Instead, to assess the accuracy with 
which these expressions predict the satellite's behavior, simulated trajectories with typical 
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Syncom n orbit elements have been run on Goddard’s particle program ’’ITEM”. These tra- 
jectories (Appendix A) confirm the validity of the derived drift equations to an accuracy well 
within the "noise levels" in the orbital elements reported for Syncom n. The equations are es- 
sentially the same as those which P. Musen has derived by a more general but more complex 
"energy perturbation" method (Reference 7). 

In Figure 1, F is a small earth- gravity perturbation force acting tangentially to an initial cir- 
cular 24-hour satellite orbit; ds is a small arc length of the satellite's path around the earth. 

At the beginning of the dynamics, the total 
energy (the sum of potential and kinetic ) of the 
satellite in a spherical earth-gravity field (Re- 
ference 8) is 

E = - 5T (1) 

where is the earth's gaussian gravity con- 
stant (3.986 x 10 s km 3 /sec 2 ). The energy added 
to the satellite by F per day is 

AE = ^Fds = 27ra s F , (2) 

where F = (1/2 rr) j) Fd» . F , in units of force 
per unit mass, is the orbit averaged energy 
perturbing force. If the orbit is purely circu- 
lar, only a tangential perturbation force can 
cause a change in the total energy. The ITEM 

simulated trajectories in Appendix A and the real Syncom n orbits both maintain eccentricities 
of the order of 0.0001 for periods up to 100 days. Equation (2) assumes the eccentricity is zero 
for the 24-hour satellite of semimajor axis a s . 

From (1), the change in energy of a 24-hour satellite is accompanied by a change in semi- 
major axis expressed by 



Figure 1 — Orbit plane of a 24-hour satellite, 
looking southerly. 


AE 


E 3 s 

2 K ) 2 ’ 


or 


Aa 


2 ( a s ) 2 AE 




( 3 ) 
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Substituting (2) into (3), the change in semimajor axis of the 24-hour near-circular orbit, per day, 
is approximately given by 


4 77 


Aa 


(*s ) 3 


(4) 


From Kepler's third law, the period of a 24-hour orbit as a function of its semimajor axis is 


H*.) 3/a 

Ts " K) 1/2 ' < 5 > 


Thus, if the semimajor axis changes by Aa s , the period change is given by 


AT 

S 


3w(a s y /2 Aa s 

Kr 


( 6 ) 


Substituting (4) into (6), the change in period, per day, of a 24-hour circular orbit is given by 


AT 


127t 2 (a s ) 7/2 F 


(7) 


The apparent net longitudinal drift rate of the 24-hour satellite's ground track with respect to the 
surface of the earth after the first sidereal day is 


(AT 1 2tt 

k (t = 1 sidereal day) = (radians/sidereal day) (8) 

s 

(See Reference 9 for example). The minus sign is taken in (8) because a gain in period is accom- 
panied by a decrease in net geographic longitude for the initially 24-hour satellite (for example, 
for the daily geographic position of the ascending node). Combining Equations (7) and (5) in (8) 
gives 


k (t = 1 sid. day) = - 12 7 (rad./sid. day) . (9) 

M a «) 

As the gain in semimajor axis is small over one day (and, in fact, small compared to a s , for 
the entire libration period of the satellite in the triaxial earth field), the drift rate will continue 
to build up linearly with time initially, adding increments of (9) each day. Thus, the net longitudi- 
nal drift acceleration of an initially 24 -hour satellite is 


k = ~ f (rad./sid. day 2 ) . 

M a s) 


( 10 ) 
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Rewriting (4) as 


4m (a s ) 3 F 

, (length units/sid. day) 


( 11 ) 


gives the expansion rate of the initially 24- hour near- circular satellite orbit due to a small but 
persistently acting orbit- averaged tangential perturbing force F . 


2. EVALUATION OF THE PERTURBING FORCE 


Figure 2 shows the position of the 24-hour satellite with respect to the earth and the celestial 
sphere. F r , F 0 and F x , earth-gravity perturbing forces in the radial, latitude, and longitude di- 
rections, are assumed to be acting on the satellite at s. Considering only the earth-gravity per- 
turbation forces arising from the ellipticity of the earth’s equator (Reference 10), Appendix B 
gives these forces as 


F 

r 


( R q/ a s ) 




cos 2 4 > cos 2 


(\ ^22)} 


( 12 ) 


NORTH 



Figure 2 — Position of a 24-hour satellite with a near-circular orbit with respect 
to the earth and the celestial sphere. 
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(^o/ 3 s) 


(12A) 


K ) 2 


|6J 22 cos 4 > sin 4> cos 2 (A. - A. 22 )} 


; (y°s ) 2 

N 2 


|6J 22 cos sin 2(^ - \ 22 )| • 


(13) 


As long as the orbit is nearly circular, F r will have negligible contribution to F. The contribution 
to F from F^ is: 


F ( ,j :i ^ cos a - K sin 1 p cos <p cos a cos 2 (X - K 22) 

K is a constant for a single orbit. 


(14) 


K 


6J 22 



In the right spherical triangle AN, s, L, note the following trigonometric relations: 


(14A) 


cos ( i ) 

tan AL 
tan 9 

(15A) 

cos a - 

tan </> 
tan 6 

(15B) 

sin 0 = 

s in ( i ) sin 9 , 

(15C) 

sin a = 

cos ( i ) 
cos 

(15D) 


From (15A), 

AL = tan -1 [tan 0 cos (i)] (16) 

Let the geographic longitude of the satellite at the ascending node (AN) be a. 0 . Counting time 
from this orbital position, the geographic longitude of the 24-hour satellite at s in its near- 
circular orbit is 


^ - X 0 + AL - &> et X Figure 2) 


or, using (16), 


A “ A Q + tan 1 [tan 9 cos ( i ) ] - 


(17) 
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w e is the earth's sidereal rotation rate. For the 24-hour satellite (starting the dynamics with 
S. at A.N., for convenience), 6 - w e t , so that (17) becomes 

A. = X„ + tan” 1 [tan 8 cos ( i )] - 8 , (18) 

approximately. The function tan" 1 [tan 6 cos (i)] -8 is even about 6 = 0 and 6 =7r/2,witha 
period of n , and behaves like a somewhat distorted sine function (Figure 3 and Appendix C). Call 
this function ax and note that for i < 33 degrees, AX is always less than 5°. Thus, using (18) and 
assuming i is sufficiently small (i < 45° proves to be a sufficient restriction on the inclination), 
cos 2 (x - x 22 ) for the 24-hour satellite can be approximated by 


cos 2 [X - X 22 ) - cos 2 (x 0 - X 22 j - 2A X sin 2 [X Q - X 22 ) cos 2y g + 2AX sin 2y 0 . (18A) 

Similarly, sin 2 (x - x 22 ) can be approximated by 

sin2(X 0 -X 22 ) + 2AXcos2(X 0 -X 22 ) = - sin 2y 0 - 2AX cos 2y 0 = sin2(x-X 22 ) (18B) 

In Figure 3, y 0 is the geographic longitude of the node of the 24-hour nearly circular satellite 
orbit with respect to the minor equatorial axis. With these expansions [(18A) and (18B)], and using 
(15B), (14) becomes 




sin d> cos 4> tan 4> r i 

K tanl? cos + ^ sin } 


Using (15C) in the above expression, the con- 
tribution to the perturbing force F due to F^ 
becomes 


sin 2 ( i ) sin 28 , 

F ( ^ - K 2 {" cos 

+ 2AX sin 2y 0 } * 

Writing ax £ ax (max) sin 28 , (19) becomes 

, r t sin 28 

F ( <£) - Ksin 2 (i) [-cos 2y Q ] 

- Ksin 2 (i)AX (max) [sin 2 y 0 sin 2 2 &] 


(19) 


( 20 ) 


Averaging F ^ over 0 < 6 <2tt, (20) gives 

tt 

sin^ iA k (max) r 

= 2 n R sin2 y 0 ] (21) 

o 


NORTH 



Figure 3— Geographic subsatellite track of 24-hour 
satellite in a near-circular orbit. 
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The contribution to F from F x is 


Fy^j n F^cos(90-a) - sin a - K sin a cos 4> sin 2y 0 - 2AX cos 2y 0 } > (22) 

from (13) and (18B). Using (15D) in (22), and noting that Ax - Ax (max) sin 2 0, as before, gives the 
contribution to F from F x as 

F (J0 = K cos ( i ) {- s in 2y 0 - 2AX (max) sin 26 cos 2y 0 } • (23) 

Averaging F (J0 over 0 < 6 <2 n, (23) gives 

F (J0 = - Kcos (i) sin 2y„ . (24) 

Thus, combining the contributions of the latitude and longitude perturbations to the average 
perturbation force over a single 24-hour orbit, (21) and (24) sum to produce 

- — — f s in 2 ( i ) AX (max) 1 

F (total) - F ( 0) + F (a.) ' - Ksin 2y 0 |cos (i) + 2 J ' (25) 

3. COMPLETION OF THE DERIVATION OF THE DRIFT EQUATIONS 
Appendix C shows that 


AX (max) - tan 1 [sec (i)] “ 45° . 

It is also shown there that, to a high degree of accuracy for i < 50°, 


sin 2 AX ( max ) 1 + cos 2 ( i ) 

cos (i) + 2 ~ — 2 ‘ 

Numerically integrated orbits have shown that the drift theory for a 24- hour satellite stemming 
from (25) is in error by more than 2 percent for i > 45°. With this restriction on orbit inclina- 
tion, using the above approximation for the inclination factor, we can rewrite the longitude drift 
and orbit expansion equations (10) and (11), evaluating F by (25), giving 


12rr 2 K [cos 2 (i) + ll „ , / , 2 \ 

y = — j — (7 2 sin 2r 0 > (rad./sid. day ) 

^E/[ a s) L J 


(26) 


( a s) 3K 




cos 2 ( 


— - sin2y 0 , (length units/sid. day) 


(27) 
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Substituting (14A) into (26) and (27) reduces these expressions to 


y ° 72 ” 2 J22 (V a .) 2 [ C ° S ' l 1 + *] sin ^0 , (rad./sid. day 2 ) 

a - - 24tt J 22 (R 0 /a s ) R 0 |^5_J Li) — -J s in 2y 0 , (length units/sid. day). 

Define a nondimensional change of semimajor axis from a s during the drift as 


Aa 

— ; so that, a, 


a 

a 


With (29A), (29) becomes 


24^ J 22 (R 0/ /a s y {— ^2 J sin 2y 0 , ( 1 /sid. day). 


Define: 


( 28 ) 

(29) 


(29A) 


(30) 


A 22 - - 72t7 2 j 22 (R 0 /aJ 2 {— S ( 2 l) + 1 j , (rad./sid. day 2 ). (30A) 

With (30A), (28) and (30) become; 

7 + A 22 sin 2y 0 = 0 , (rad./sid. day 2 ) (31) 

A 22 sin 2y 0 

a 1 ~ 3^7 = 0 , (1/sid. day). (32) 

Note that n in (32) has dimensions of rad./sid. day. It must be understood that (31) describes 
the net daily geographic acceleration of the initially 24- hour satellite with respect to the earth's 
minor equatorial axis. Stated another way, (31) describes the geographic drift of the entire 
originally stationary, figure - 8 ground track (Figure 3). Similarly, (32) describes the net daily 
orbit- expansion rate of the 24-hour satellite. In particular, it is convenient to treat the motion of 
the ascending node of the orbit in geographic longitude as a reference for the entire configuration. 
In what follows, therefore, 7 will refer always to the geographic longitude of the ascending node 
east of the equatorial minor axis; > 0 will refer to the initial geographic longitude of the A.N. east 
of the minor axis, at the start of the dynamics under consideration. (31) and (32) can thus be re- 
written in terms of the general nodal longitude position y, to give the relevant partially uncoupled 
long-term drift and orbit- expansion differential equations for the near-24-hour near-circular 
orbit satellite: 


7 + A 22 sin 27 = 0 , (rad./sid. day 2 ) 


(33) 
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A 22 sin 2 y 
dm 


0 , (1/sid. day). 


(34) 


4. GENERAL CONSIDERATIONS OF THE SOLUTIONS OF THE DRIFT EQUATIONS 
Equation (33) can be integrated directly for the geographic drift rate by noting that 

■■ _ dy d(>) 2 d(y) 2 

7 " dt 2y dt W ' 

Thus (33) can be separated to 


d(y) 2 - - 2A 22 sin 2y dy . (35) 

Since the variables (y ) 2 and y are separated in (35), (35) integrates to 

(r ) 2 = A 22 cos 2y + c x . (36) 

With the initial condition that 7 y 0 at y = y 0 , (36) becomes 

y = [(^o) 2+A 22 (cos 2y - cos 2r 0 ) J 1/2 , (37) 

giving the drift rate of the 24 -hour satellite as a function of the initial drift rate y 0 , the earth- 
gravity constant A 2 2 , the initial longitude east of the minor axis y 0 , and the instantaneous longi- 
tude 7. Returning to the semicoupled system of equations (33) and (34), the explicit dependence of 
the equations on the location from the minor axis and the magnitude of the equatorial ellipticity 
may be eliminated by multiplying (33) by 1/3^ and adding the resulting equation to (34), giving 

y + 3tt a, = 0 , (38) 


(38) can be rewritten as 

"dT^ + 3-77 3j = 0 - d(y) + 3tt a j dt = d(y) + 37rd(aj) . (39) 

Separation of the variables y and a t is thus achieved in (39). (39) integrates directly to 


377 3J + y - C 2 . (40) 

With the initial conditions: a 1 = 0, when 7 = 0 (the satellite is in the momentarily stationary 
ground-track configuration); C 2 = 0. If y 0 in (37) is also the longitude of this initially stationary 
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orbit, (y 0 ) 2 = 0 there, and (37) in (40) yields for a t , the semimajor axis change from "synchron- 
ism" in the drift motion, 


( cos “ cos 2y 0 ) w 
3m 


(41) 


(41) shows explicitly that the semimajor axis is bounded in long-term drift from a stationary 
orbit. From (33), since A 22 > 0, if 0 < y 0 < 90°, drift proceeds towards the nearest longitude of 
the earth's equatorial minor axis (in a - y direction). If -90°< y 0 < 0°, (33) shows that drift again 
proceeds toward the nearest minor axis longitude (in a +7 direction). Thus, in all cases of drift 
from a stationary geographic configuration, cos 2 y - cos 2 y 0 is a positive function which has a max- 
imum when 7=0 (when the satellite has drifted over the longitude of the minor axis). Thus (41) 
gives (for the librations of a 24-hour satellite) 


Sj (max) 


(A 22 ) 1/2 (l~ cos 27 0 ) 1/2 

3tt 


(42) 


Again it is noted that -n in (42) has units of rad/sid. day. An absolute maximum semimajor axis 
change in the drift occurs when the "synchronous" condition is established near the longitude of 
the major equatorial axis. Here, y 0 = -90°, cos 2y 0 = -1, and 


(2A 22 ) 1/2 

aj (absolute maximum for a librating 24- hour satellite) = — 3 - — (43) 

For the constant J 22 = - 1.7 x 10 " 6 (derived in this study from long-term observations on the 
drift of the Syncom n satellite) and using the additional constants from this study: i = 33°, 
a s - 42166 km, R 0 = 6368.388 km; (30A) gives 

a 22 = 23.2 x 10 ~ 6 (rad./sid. day 2 ) . 


(43) then gives 


a j (absolute max.) = .72 x 1(T 3 , from which, by (29A), 

Aa (absolute max. from a "synchronous" condition near 

the equatorial major axis, for a satellite of i = 33°) = 30.7 km . 


Thus the assumption made in (10) and (11), to approximate the slightly varying semimajor axis by 
a s (a constant) throughout the drift motion, appears amply justified. 

Figure 4 is a graph of (41) for a, vs. 7 (the longitude with respect to the nearest minor axis 
location) as a function of y 0 , the longitude in the initially stationary configuration. 
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Figure 4— Libration with longitude of the semimajor axis of a 24-hour satellite as a 
function of the longitude of the initially stationary configuration. 


Note that (41) allows equal ± values for for each > . Suppose the satellite is initially at 'y 0 
(position 1 in Figure 4) from the nearest location of the minor axis: From (33), sin 2y 0 being posi- 
tive, the satellite begins to drift west (attaining a negative drift rate) towards the minor axis. But, 
from (40), since C 2 = 0, a 1 = - (y)/37i> 0; the drift therefore proceeds counterclockwise in Fig- 
ure 4, around the central point of the minor axis and a 1 = 0 , along the upper portion of the two- 
valued arc determined from (41). 

The same situation holds for the motion beginning or stemming away from the "synchronous" 
longitude at ~y 0 , position 2 in Figure 4. Here sin 2y 0 is negative, and the drift proceeds at a posi- 
tive rate to the east. Again from (40), as soon as the satellite leaves position 2, a x - - (y)/3n < 0, 
and the circulation continues in a counterclockwise direction. Every trajectory in the phase plane 
a^y may be conveniently defined by the constant C, of the "energy integral" of the drift motion 
(36). Since (33) is the equation- of- motion defining the large-angle oscillations of a mathematical 
pendulum (in the case of the 24-hour-orbit satellite, the point of symmetry is the minor axis where 
2y = 0), it can be expected that the general solutions in that theory apply to the long-term libra- 
tions of the "synchronous satellite” (Appendix D). For example, in (36), with a momentarily 
"synchronous" condition at y 0 being given by y 0 = 0, the "energy constant" is evaluated as 


Cj ^22 2 y q . 
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With this evaluation, (36) becomes 


(r) 2 “ A 22 cos 2y ° - A 22 cos 2y 0 . 


(44) 


Solving for the initially "synchronous" longitude as a function of any longitude in the drift and the 
corresponding longitude rate, (44) gives 


y 0 


cos' 1 I cos 2 y - 


(r) 2 ~ . 
a 22 ^ 


(45) 


Since (y)J/. A 22 > 0, the argument of cos -1 in (45) is always less than or equal to 1. Thus, as long 
as cos 2 y - (y)//A 32 > - i, (45) will give a real solution for the momentarily "synchronous" longitude 
with respect to the minor axis. But, if cos 2 y - (y) 2 / A 22 < - l, there will be no real momentarily 
"synchronous" configuration for the near-24-hour satellite. With this energy, the world- 
circulation regime commences, corresponding to the over-the-top, high-energy regime of the 
mathematical pendulum (Reference 3). The above inequality implies that, for the commencement 
of "world circulation" for the near-24-hour satellite, 


(r) 2 i A 22 (1 + cos 2y), or 

(r) 2 > 2A 22 cos 2 y . (46) 

When 2 y = 0, or the satellite is over the minor axis, (46) allows the maximum possible drift rate 
for a librating 24-hour satellite: 


•y( max ) fo , juration “ ( 2A 22 ) V2 , (rad./sid. day) . (47) 

For example, using the reported value of A 22 = 23.2 x 10‘ 6 rad./sid. day 2 for the inclination of 
the Syncom II satellite, (47) gives 


y(max) 

for libration with J 22 


-1.7 X 10' 6 , 


33<) = (46.4 x 10~ 6 ) v2 = .39 degrees/day. (48) 


5. APPROXIMATIONS TO THE EXACT DRIFT SOLUTIONS FOR PERIODS VERY CLOSE TO 
SYNCHRONOUS 

Expanding the drift from the "synchronous" longitude (y - y g in this section) in a Taylor 

series, with respect to increments of time At from the momentarily stationary condition, 


X At ) 


>0 


To At 


+ .. ( A t ) 2 

+ y 0 2 


+ Z'o 


( A t ) 3 . [ 4 ] (At ) 4 

““ 6“ T >o 24 


[ 5 ] ( At ) 


m~ + y 0 


[ 6 ] 


(At ) 6 

~7W 


•M 


(At ) 7 

^0 




yo 40320 


(50) 
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Differentiating (33) six times with respect to time, it is clear that all derivatives in (50) can be 
written as functions of A 22 , y 0 and y 0 . Noting that y(At) - y 0 - A\ (the geographic longitude with 
respect to the "synchronous" configuration) and y 0 = 0, (50) can be shown to reduce to the 
expansion 


A\ = (-A 22 sin 2y 0 ) ^- + [(A 22 ) 2 sin4y 0 ] ^- + [(A 22 ) 3 sin2y 0 (4sin 2 2y 0 - l)] 


- [(A 22 ) 4 sin4y 0 (34 sin 2 2y 0 - l)] 


8 


(At) 

10080 


(51) 


It is apparent that, as At - 0, the higher order terms of (51) become increasingly more insignifi- 
cant to the total drift, in comparison to the terms of lower order. 

In Appendix D, the exact "elliptic integral" of motion from (33) is calculated from a syn- 
chronous longitude of 60° east of the minor axis. This calculation demonstrates that the simple 
term- inclusion-time criterion below gives an adequately converging series to the "exact" drift. 

In the actual reduction, all higher order terms in (51) which are less in magnitude than the root 
mean square (rms) error of the observed Syncom n longitudes, are ignored. Section 7 of this re- 
port shows that this rms error of longitude determination for the ascending equator crossings of 
Syncom n from August 1963 to March 1964 has been of the order of ±.025 degrees. Thus, 0.025° 
is used below in forming the minimum-time-term-inclusion criterion for each term of (51). 

A 22 is assumed to be 23.2 x 10~ 6 rad./sid. day 2 . 

A) For inclusion of the (At) 4 term: 


sin 4y 0 is maximum when y 0 = ±22.5° and ±67.5° . 


\2 (At)’ . 
( A 22j 24 


Therefore, 

|AX majt (from the fourth-order term)| 

Solving (52A) for At, when |AA. max (4th order)| = .025°, 

At (min. fourth-order term inclusion) = ^.025 x 24/57.3 X [23.2 x 1(T 6 ] 2 ) 


(52A) 


1/4 


= 66.5 sid. days from "synchronism" . 


B) For inclusion of the (At ) 6 term: 


sin 2 y 0 (4 sin 2 2y 0 - 1 ) | is maximum when y 0 = ±45° . 


Therefore, 


A^max (from the sixth-order term)| 


( A 22 ) 3 


(At) 6 

' 


(52B) 
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Solving (52B) for At , when (Aa.^ (sixth order) j = .025°, 

At (min. for sixth-order term inclusion) = (.025 x 60/57.3 x [23.2 x 10~ 6 ] 3 ) 

= 113. sid. days from "synchronism" . 

C) For inclusion of the (At) 8 term: 

| sin 4 y g (34 sin 2 2y 0 - l) | is maximum when y 0 = ±59 V4° . 

Therefore, 

j AA majt (from the eighth-order term)] = 21.2 (a 22 )‘ i mogp 
Solving (52C) for At, when |AA max (eighth order) | = .025°, 

At (min. for eighth-order term inclusion) = (l0080 x .025/21.2 x 57.3 x [23.2 x lO -6 ] 4 ^ 

= 171. sid. days from "synchronism" . 

Similarly, expanding a, (t) in a Taylor series about the time of "synchronism", (t 0 , y 0 ): 
(*o + At ) = ( a i ) 0 + (a^oCAt) + (a,)o ^ + ( 3 .)o ^ ♦ •" ' 

But, from (34), 


Aj 2 sin 2 y 0 
( a l)o " 377 


Differentiating (5 3 A) with respect to time, 

... . _ ^22 COS 0 
( 3 l) 0 377 


since y 0 = 0. Differentiating (53B) with respect to time, 

... _ ~ 4 (/q) 2 A 22 sin 2 >o 2 Jq A 22 cos 2y 0 _ - (A 22 ) 2 sin4y 0 

( a i/0 377 + 377 ~ 377 


using equation (33). From the conventional definition of a,, (a,) 0 = o. (53) then becomes 

/A . , (A 22 sin 2y 0 ) At (A 22 ) 2 sin4 7o (At) 3 

a t (at At from "synchronism") = + ... , 

with the results of (53 A), (53B), and (54) in (53). 


(52C) 


(53) 


(53A) 


(53B) 


(54) 


(54A) 
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Section 7 shows that the rms error of semimajor axis determination for Syncom II (including 
sun and moon "noise") is of the order of ±0.5 km. Therefore, the rms error to be expected in a t 
is of the order of .5/42166. = 1.185 x 10' 5 . Following the procedure for the longitude drift, 

1.185 x 10' 5 is used below to determine the minimum time for the inclusion of the terms beyond 
the first on the righthand side of (54A), to ensure adequate convergence of the infinite series for 
3j (At). 

A) For inclusion of the (At) 3 term: 

| sin 4y 0 | is maximum when y 0 = ±22.5° and ±67.5° . 


Therefore, 

j a i [from the third-order term of (54A)]| = 
Solving (55) for At, when |aj (max)j = 1.185 x 10' s , 

At [min. for the third-order term inclusion in (54A)] = 


From a "synchronous" configuration at 54.8° west of Greenwich, on or about 6 September 1963, 
Syncom n drifted to 59.2° west of Greenwich on 28 November 1963, where it was "stopped" by the 
tangential firing of on-board cold-gas jets. A second free-drift period followed from a "synchron- 
ous" configuration at 59.2° west on about 29 November 1963, to 66.3° west on 18 March 1964, 
where the on-board tangential jets were fired to speed up the westward drift. Of the 34 separate 
orbits calculated by the Goddard Data and Tracking Systems Directorate for these free-drift 
periods, only 7 fell outside the minimum 66-day period around a condition of "synchronism", for 
which the inclusion of higher order terms in (54A) would be necessary in reducing the drift data 
according to that theory. The data reduction of Section 7 includes only those orbits falling within 
the minimum 66-day period around "synchronism". Further refinement of this reduction to in- 
clude the 7 outside- of- synchronous orbits (according to the criterion of this chapter), will be made 
in the near future. This refinement is not expected to materially affect the results of this report. 


(A 22) 2 (At ) 3 


18 /f 


(55) 


1.185 X 10' 5 X I877 


23.2 x 10' 


1/3 


108 sid. days from 
synchronism . 


6. DETERMINATION OF EARTH- EQUATORIAL ELLIPTICITY FROM TWO OBSERVATIONS 
OF DRIFT ACCELERATION AT A GIVEN LONGITUDE SEPARATION 

Given two independent near- synchronous drifts (in the sense discussed previously), whose 
momentarily synchronous longitudes (y 0 ) 1 and (y 0 ) 2 are separated by v\. Let the two drift 
accelerations at these two "synchronous" configurations be (y 0 ) t and [y 0 ) 2 • The drift ac- 
celerations may be determined from drift-data reduction according to the theory of (51). 
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From (33), 


(r'o)i " - (A 22 )j sin 2(y 0 )i 


(56) 


(^0)2 (^22)2 s * n ^[(^0)1 • 

since ( y 0 ) 2 - (y,,^ - — - = (\ 0 ) 2 - {^ 0 ) 1 . Expanding (57) and dividing by (56) gives 


(^0)2 

( A 22 ) 1 

(^’0)1 

(^22) 2 


cos 2VA + sin 2V\ cot 2^y 0 )j . 


Solving (58) for (y 0 ) j , 


(y 0 ), = 2 tan_1 i 


sin (VA-) 


(■y'o) 2 (^22) 1 
(t'o) 1 (^22)2 




- cos 2 (V\) 


(57) 


(58) 


(59) 


The quadrant of (y 0 ) 2 is either the first or the fourth, because drift acceleration is always in 
the direction of the nearest longitude extension of the earth’s minor equatorial axis. Once the 
minor axis is located by (59), the absolute value of J 22 in the earth's triaxial gravity field can be 
determined through (56) and (30A), for example, as 


( A 22j 

li 

(y'o)i 




72tt 2 | 

[Vl*.).] 

^ l 

cos 2 (i)j + 1 

72n 2 sin2(y 0 )j [R 0 /( a s)i] 

1 2 

' j 

cos 2 ( i ) t + ll 

(60) 

L 2 

2 J 


Note that the units of (y 0 ) 2 in (60) must be those of radians/sidereal day 2 so that J 22 will be di- 
mensionless. Note also that in (59), using the result of (30A), 


( A 22) 1 

( 3 s) 2 

2 

cos 2 ( i ) j + 1 

( A 22 ) 2 



cos 2 ( i ) 2 + 1 


(61) 


Using (59), since (<v 0 ) 2 is known from the data reduction (the geographic longitude of the "syn- 
chronous" configuration), the geographic longitude of the nearest minor axis location can be 
calculated as 

^22 ' (^0)1 ~ (^0)1 • ( 61 A) 

Similarly, the geographic longitude of the nearest major equatorial axis location can be calculated 
from 

^22 = (\>)i “ Wi + 90 ° • (6iB) 


(See figure 3). 
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Following the theory of Reference 10, the difference in major and minor equatorial radii of 
the earth's triaxial geoid (a 0 - b 0 ), is related to the gravity constant J 22 by 

a 0 ' b 0 ~ _ 6R 0 J 22 • (62) 


7. REDUCTION OF 27 SYNCOM II ORBITS TO DETERMINE THE EARTH'S EQUATORIAL 
ELLIPTICITY 


Appendix A tabulates the 27 Syncom n orbits from which the reduction below was made. 

Table 1 gives the estimated ascending equator crossings nearest to the epoch of these orbits. 

These were calculated by hand, and therefore are listed only to 0.01 degrees and 0.01 days. The 
technique used was to locate from the Nautical Almanac, the geographic longitude of the ascending 
node at epoch through the reported right ascension of the ascending node for the orbit, and the 
hour-angle of the vernal equinox calculated at epoch. The geographic longitude of the ascending 
equator crossing was then estimated by turning the earth back through the orbit angle from the 
ascending node to the satellite at epoch. This latter quantity was estimated as w - M for the near- 
circular orbit of Syncom II. A correction factor to this orbit angle — the ratio of the satellite's 
period to the earth's sidereal period — was applied for orbits whose period was sufficiently dif- 
ferent from the earth's. The nodal longitude at epoch, minus this reduced nodal excursion angle, 
is the estimated "ascending equator crossing nearest to epoch" reported in Table 1. (See Appendix 
E for an example of this calculation.) 

Table 2 gives the Goddard- reported semimajor axes for these 27 orbits. Truncating equa- 
tions (51) and (54A) at their first righthand terms: 


AX. (longitude drift from "synchronism”) 


(a 22 


sin 2 y 0 ) 


(At) 2 


(63) 


sin2y 0 (At) 

a t (semimajor axis change from "synchronism”) a A 22 ^ • (64) 


Let the drift time be given from a certain arbitrary base time by T. Let T 0 be the time of 
"synchronism" from the base time. Let the drift be given from a certain arbitrary geographic 
longitude by X . Let x q be the geographic longitude from this base longitude, of the "synchronous" 
configuration. Then: 


At = T - T 0 , and 


AX = x - x 0 


; ) : 


With these changes, (63) and (64) become (noting that a t = 

_ M J Sln 2 X 0 ) 

A - A 0 - 2 l T 2TT 0 V) • or 

T 0 ( A 22 S ^ n 2 ^o)1 f 1,1 A 22 S ^ n 2 ”Xo 

2 h T 1 T o A 22 »in2y 0 K T 2 


A = i A 0 - 


(65) 
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Table 1 


Estimated Ascending Equator Crossings Nearest the Epoch of 27 Syncom II Orbits. 


First Drift 
Orbit # 

Time from 
20.0 Aug. 1963 
(days) 

Ascending Equator 
Crossing in Degrees 
West of 50.0° West 

Second Drift 
Orbit # 

Time from 
26.0 Nov. 1963 
(days) 

Ascending Equator 
Crossing in Degrees 
West of 50.0° West 

1-1 

2.12 

4.89 

2-1 

1.86 

9.17 

1-2 

7.11 

4.83 

2-2 

7.84 

9.17 

1-3 

11.09 

4.78 

2-3 

13.83 

9.22 

1-4 

16.08 

4.74 

2-4 


9.38 

1-5 

20.07 

4.77 

2-5 

41.75 

10.15 

1-6 

23.06 

4.78 

2-6 

44.74 

10.36 

1-7 

28.05 

4.85 

2-7 

55.71 

11.02 

1-8 

31.04 

4.90 

2-8 

64.69 

11.76 

1-9 

38.02 

5.06 

2-9 

71.67 

12.32 

1-10 

42.01 

5.09 

2-10 

76.66 

12.81 

1-11 

48.99 

5.45 

2-11 

83.64 

13.49 

1-12 

54.97 

5.68 




1-13 

62.95 

6.09 




1-14 

70.93 

6.60 




1-15 

77.91 

7.14 




1-16 

83.90 

7.61 





-100.0 

First free drift 






period ends at an 






ascending equator 






crossing of -9.15° 






west of 50.0° west. 





Table 2 


Goddard-Reported Semimajor Axes for 27 Syncom n Orbits. 


First Drift 
Orbit # 

Time from 
20.0 Aug. 1963 
(days) 

Semimajor Axis: 
(42160.0 + Data; km) 

Second Drift 
Orbit # 

Time from 
26.0 Nov. 1963 
(days) 

Semimajor Axis: 
(42160.0 + Data; km) 

1-1 

2.27 

4.58 

2-1 

2.04 

5.89 

1-2 

6.71 

4.52 

2-2 

8.00 

7.20 

1-3 

11.00 

6.02 

2-3 

14.00 

7.18 

1-4 

16.00 

6.39 

2-4 

20.71 

8.17 

1-5 

20.00 

6.35 

2-5 

41.71 

8.01 

1-6 

23.08 

6.55 

2-6 

44.25 

9.90 

1-7 

28.08 

6.70 

2-7 

55.88 

ll.vJ 

1-8 

31.08 

7.42 

2-8 

64.83 

11.91 

1-9 

38.08 

7.51 

2-9 

71.67 

12.89 

1-10 

42.08 

8.88 

2-10 

76.79 

13.31 

1-11 

49.08 

9.14 

2-11 

83.71 

14.89 

1-12 

55.08 

9.78 




1-13 

63.08 

11.51 




1-14 

71.00 

11.09 




1-15 

7 8.' 00 

12.15 




1-16 

84.21 

12.51 





-100.0 

First free drift 
period ends with a 
semimajor axis of 
-42174.5 km. 
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(Note that from (65), X = - A 22 sin 2y 0 = % from (33). This result is valid only for orbits sufficiently 
close to "synchronous,” as discussed previously.) 


a 


+ a 


A 22 sin2y 0 


377 


[t - T fl 


, or 


T 0 A 22 sin2 >o\ (a s A 22 sin 
a - a. | 1 ■ )+T 


3tt 

(66) and (65) may be written with determinable coefficients as 

A. = d Q + dj T + d 2 T 2 , 

f 

a " e 0 + e, T , 


377 


where: 


From (69), 


T 2 A 22 sin 2y 0 
~ — o 2 

d, = (+ A 22 T 0 sin 2y 0 ) , 

A 22 sin 2y 0 


T 0 ^22 s * n ^Xo 
37 


sin 2y 0 

e 1 ~ a s A 22 3n 


( 66 ) 


(67) 

( 68 ) 


(69) 


T 0 = - d/2d 2 , 

N 2 


4d, 


Xo 


A 22 sin 2y 0 - 2d 2 


(70) 

(71) 

(72) 


Alternately, and as an internal check on the theory of the coupling of the drift and orbit expansion, 

~3me , 

J 0 = “ A 22 s * n 2y 0 " ■ 


— , implying 


3-ne 1 j 2a 


(73) 
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In equation (73), the units of d 2 must be radians/sidereal day 2 , and the units of e, must be length/ 
sidereal day so that the equation will be dimensionally correct. The semimajor axis at the "syn- 
chronous" configuration is calculated from (68) for T = T 0 : 

a s D e o + e if 0 • (74) 

For the first drift period (orbits 1-1 through 1-16), the best estimates (in the "least squares" 
sense) of the coefficients (d)j and (e)j, obtained by fitting (67) and (68) to the data in Tables 1 and 
2, have been found to be: 

(d 0 )i = (4.941 ± .018) degrees 
( d i)i = -(.0216 ± .0010) degrees/solar day 
(d 2 )j = (6.37 ± .11) x 10 -4 degrees/solar day 2 
= (6.33 ± .11) x 10~ 4 degrees/sid. day 2 

( e o)i = (4-35 ± .19) km 
( e i)i = (-0993 ± .0^42) km/solar day 
(.0990 ± .0042) km/sid. day . 

The mean value of the inclination during this period was 

(i), = 33.018 ± .005 degrees . 


From (70), 

(To)! = ^16.95 ^l;Sfc)days from 20.0 August 1963 . 

From (71), 

( A c )i = (4.76 ± .03) degrees west of 50.0 degrees west longitude . 

From (72), 

(y 0 )! = -(1.27 ± .02) x 10" 3 degrees/solar day 2 = -(2.20 ± .04) x 10" 6 rad./sid. day 2 . 
From (74) and the above value of (t q )j , 


(ajj = (42166.0 ± .2) km . 
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For the second drift period (orbits 2-1 through 2-11), the best estimates (in the "least squares" 
sense) of the coefficients (d) 2 and (e) 2 , obtained by fitting (67) and (68) to the data in Tables 1 and 
2, have been found to be: 

(d 0 ) 2 = (9.156 ± .017) degrees 
(dj) 2 = -(.0030 ± .0010) degrees/solar day 
(d 2 ) 2 = (6.59 ± .11) x 10' 4 degrees/solar day 2 
= (6.55 ± .11) x 10~ 4 degrees/sid. day 2 

(e 0 ) 2 = (5.70 ± .42) km 
(ej) 2 = (.0994 ± .0080) km/solar day 
= (.0990 ± .0080) km/sid. day . 


The mean value of the inclination during this period was 

( i ) 2 = 32.851 ± .010 degrees . 

From (70), 

(t 0 ) 2 = 2.3 ± .8 days from 26.0 November 1963 . 

From (71), 

(\ 0 ) 2 = 9.15 ± .02 degrees west of 50.0 degrees west longitude . 

From (72), 

(-y 0 ) 2 = -(1.32 ± .02) x 10' 3 degrees/solar day 2 = -(2.29 ± .04) x 10 -6 rad./sid. day 2 . 

From (74) and the above value of (t 0 ) 2 , 

(a s ) 2 = (42165.9 ± .4) km . 

(See Figure 5 for a graph of this orbit data and reduction for the two drift periods.) Combining 
the above results of the two free-drift periods, from (61), 


(^22) 1 


(a 22 ) 2 


(42165.9 ± .4/42166.0 ± .2 ) 2 
.99845 ± .00014 . 


[cos 2 (33.018 ± .005)+ l] 
[cos 2 (32.851 ± .010)+ 1 ] 
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The longitude separation between the two drift periods is given by 


VV = (\>) 2 - (^ 0)1 = [-(59.15 ± .02)] - [-(54.76 ± .03)] degrees 

= -(4.39 ± .05) degrees geographic longitude . 


Thus 


UK = -(8.78 ± .10) degrees geographic longitude . 


Therefore, from (59), the location of the minor equatorial axis with respect to the "synchronous" 
longitude during the first free-drift period (54.76 ± .03 degrees west of Greenwich) is 


Mi 


o tan 


sin [-(8.78 ± .10) 

(1.32 ± .02) ( j99845 ± 00014) ~ cos 
(1.27 ± .02) 



[-(8.78 ± .10)] 



degrees east of the minor equatorial axis . 


From (61B), the best estimate of the location of the major equatorial axis is 

/ + 4 \ / + 4 \ 

\ 22 = -55 - 154 _gJ + 90 = -119 _gl degrees geographic longitude . 

From (60), the best estimate of the triaxial gravity coefficient J 22 is 


-(2.20 ± .04) x 10' 6 

J 22 “ 

12n 2 

sin 2^54 + jL) 

(6378.2/42166.0 ± ,3) 2 

cos 2 (33.018 ± .005) + 1 

2 


= -(l.67 +-JJJ) x 10‘ 6 . 


The mean equatorial radius, taken as R 0 = 6378.2 km, is a compromise between a number of cur- 
rently used values. It is stated above without error. The likely error in (a s ), has been increased 
arbitrarily by 0.1 km. to account for the likely uncertainty in R 0 . 

Using the above estimate of J 22 from observations on Syncom II drift, the difference between 
the major and minor equatorial radii of the triaxial geoid is, by (62), 

a 0 - b 0 = 64 ^ ^ meters = 210 ^^feet . 

Comparing the deviation due to earth ellipticity with other higher order earth- gravity devia- 
tions (Appendix B and Reference 10), we note that the above figure implies a maximum deviation 
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from the mean earth sphere, due to the ellipticity of the equator, of 


AR 0 = 105 feet • 
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LIST OF SYMBOLS 


J 


nm’ 


\ 


nm 


F 

e 


a, a 

S 


ds 




T , T 

p 


S 


K, r, <p 


(').<“).() 


R 


0 

1 




0 


A ( ) 
^0 


y 


a 


l 


Spherical harmonic constants (order n, power m) of the earth's gravity potential 
A gravity force per unit mass acting on a 24- hour satellite 

(except in Appendix F) The argument from the ascending node to the satellite position 
for the 24- hour orbit 

The instantaneous semimajor axis, and the "momentarily synchronous" semimajor axis 
of the orbit of the 24-hour earth satellite. (a s , estimated to within 2 km, is 42166 km.) 

A small arc length of a space trajectory 

The earth's gaussian gravity constant (3.986 x 10 5 km 3 /sec. 2 ) 

The orbital period for a satellite, and the "momentarily synchronous" period of a 24- 
hour satellite (i.e., the earth's sidereal rotation period) 

Geographic longitude, geocentric radius, and geocentric latitude of the 24-hour satellites 
position 


^7=4 , d ^ ) , - d ^ ^ : time differentials 
* dt 2 dt" 

The mean equatorial radius of the earth (6378.2 km) 
The inclination of the orbit of the 24-hour satellite 


The "initial" geographic longitude of the satellite, or the ascending node of the 24-hour 
satellite's orbit at the start of the dynamics under consideration 

The earth's sidereal rotation rate (.7292115 x 1CT 4 rad./sec.) 

Real time 

A small argument ( ) 

The geographic longitude (positive to the east) of the ascending node of the 24-hour satel- 
lite's orbit with respect to the earth's minor equatorial axis' longitude location, at the 
start of the dynamics under consideration 

The geographic longitude (positive to the east) of the 24- hour satellite, or the ascending 
node of the satellite's orbit with respect to the longitude of the earth's minor equatorial 
axis 

(a -a s )/a s ; a nondimens ional semimajor axis change for the 24-hour satellite's orbit, 
with respect to the "momentarily synchronous" semimajor axis 
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A 22 The driving function causing drift and orbit expansion of a 24-hour satellite in a "tri- 
axial" earth-gravity field; a constant for a given 24-hour orbit inclination 

( ) The argument ( ) at the start of the dynamics under consideration 

( ) The argument ( ) at a specified location n (except in Appendix A: ( )s; the argument 

for the simulated trajectory) 


vx The geographic longitude difference between two "momentarily synchronous" 24-hour 
satellite configurations 

( a o +b o 

a 0 , b 0 The major and minor equatorial radii of the "triaxial" earth I R 0 = - — 5 — , according to 
the definition in Reference 10^ 

w The argument of perigee in a satellite orbit: The orbit angle (from the center of the 
earth) from the ascending node to perigee 


M The mean anomaly of the satellite in its orbit: The orbit angle (from the center of the 
earth) from perigee to a point M in the orbit, where M = t being the real time 
since perigee passage and T p the period of the satellite's orbit 

(e n ) , (d n ) Determinable coefficients in the drift and orbit- expansion equations (67) and (68) 

T 0 The time of "synchronism" from an arbitrary base time of reckoning T 


F( i ) The inclination factor in the triaxial driving function A 22 
V E The gravity potential of the earth 

g 0 , g s The radial acceleration of the earth's gravity field at the earth's surface, and at the 
altitude of the "synchronous” satellite 


m A test mass 

F(k, £) The elliptic integral of the first kind with argument (or amplitude) £ and modulus k 
A The longitude location of the vernal equinox 
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Appendix A 

REDUCTION OF SIMULATED PARTICLE TRAJECTORIES FOR 
EARTH EQUATORIAL ELLIPTICITY 

Tables A-l and A-2 present data taken from a numerically integrated particle trajectory of a 
triaxial earth in the presence of the sun and moon's gravity field. Only perturbed equations of 
motion from a periodically rectified Keplerian reference orbit are actually integrated by the 
digital computer program (called ITEM at the Goddard Space Flight Center). For the 3 months' 
real orbit time of these trajectories, the accumulated truncation and roundoff error in the 
numerical integration is believed to be negligible for the purposes of this reduction. The ini- 
tial position and velocity conditions for these simulated trajectories were the same as those re- 
ported for the "actual" Syncom n orbits 1-2 (for the trajectory of Table A-l) and 2-3 (for the 
trajectory of Table A-2). The program used the earth gaussian- gravity constant 

= 3.9862677 x 10 s km 3 /sec 2 , 

which is the gravity constant used by the GSFC Data and Tracking Systems Directorate in comput- 
ing the elements of satellite orbits from radar and Minitrack observations. The best estimates 
(in the "least squares" sense) of the coefficients (d) sl and (e) sl , obtained by fitting the drift and 
orbit expansion equations (67) and (68) to the data in Table A-l, have been found to be 

(d 0 ) sl = 4.841 ± .004 degrees 
(d,) sl = -(1.22 ± .03) X 10 -2 degrees/solar day 
(d 2 ) s i = (6.303 ± .038) x 10 " 4 degrees/solar day 2 
= (6.268 ± .038) x 10' 4 degrees/sid. day 2 
(e 0 ) sl = 5.45 ±41 km 
( e i) s i = (-091 ± -010) km/solar day 
= (.091 ± .010) km/sid. day . 

The mean value of the inclination during this first simulated trajectory period was 

( i ) sl = 33.005 ± .003 degrees . 
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Table A-l 


Data from Simulated Trajectory Beginning with the Elements of Syncom II Orbit 1-2. 


(J 22 = -1.68 x 10 6 ,R 0 = 6378.388 km, y 22 = -108.0: Input into Trajectory Program) 


Time from 
26.709 Aug. 1963: 
(Solar Days) 

Ascending Equator 
Crossing: 
(Degrees West of 
50.0° West Geog. Long.) 

Semimajor Axis 
(42160.0 + Data; km) 

Inclination 

(32.0 + Data; Degrees) 

2.390 

4.816 

5.27 

1.089 

8.374 

4.783 

7.09 

1.072 

14.358 

4.792 

6.01 

1.056 

20.341 

4.861 

8.12 

1.043 

26.324 

4.954 

7.13 

1.025 

32.308 

5.101 

8.98 

1.019 

38.292 

5.291 

8.31 

.997 

44.276 

5.537 

9.67 

.991 

47.268 

5.678 

10.38 

.983 

50.260 

5.821 

10.03 

.972 

53.253 

5.975 

9.42 

.967 

56.245 

6.144 

9.74 

.966 

59.237 

6.326 

11.09 

.960 

62.229 

6.522 

11.94 

.957 


Table A -2 


Data from Simulated Trajectory Beginning with the Elements of Syncom II Orbit 2-3. 


( J 22 = -1.68 x 10~ 6 , R„ = 6378.388 km, y 22 = -108.0: Input into Trajectory Program) 


Time from 
10.0 Dec. 1963: 
(days) 

Ascending Equator 
Crossing: 
(Degrees West of 
50.0° West Geog. Long.) 

Semimajor Axis 
(42160.0 + Data; km) 

Inclination 

(32.0 + Data; Degrees) 

0.823 

9.243 

6.88 

.881 

5.809 

9.351 

7.30 

.881 

10.796 

9.495 

9.41 

.877 

15.783 

9.666 

8.15 

.864 

20.769 

9.885 

9.85 

.864 

25.756 

10.134 

10.60 

.850 

30.743 

10.401 

9.95 

.842 

35.730 

10.708 

11.81 

.841 

40.717 

11.044 

12.18 

.825 

45.704 

11.412 

11.58 

.816 

50.692 

11.830 

13.93 

.808 

55.679 

12.259 

13.11 

.790 

58.672 

12.534 

13.07 

.785 

60.667 

12.724 

13.69 

.784 


From (70), 


( T o) s i = 9.68 ± .30 days from 26.709 August 1963 
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From (71), 


(-0) Si 


4.782 


+ .008 
- .007 


degrees west of 50.0 degrees west longitude . 


From (72), 


(^o)si = -(1-261 ± .008) x 10' 3 degrees/solar day 2 
= -(2.188 ± .013) x 10~ 6 rad./sid. day 2 . 

From (74) and the above value of (t 0 ) sl , 

(«,)„ = 42166.3 ± .4 km . 

The best estimates (in the "least squares" sense) of the coefficients (d ) s2 and (e ) s2 , obtained by 
fitting the drift and orbit expansion equations (67) and (68) to the data in Table A- 2, have been 
found to be 


(d 0 ) s2 = 9.224 ± .004 degrees 
(dj) s2 = (1.830 ± .028) x 10" 2 degrees/solar day 
(d 2 ) s 2 = (6.501 ± .042) x 10“ 4 degrees/solar day 2 
= (6.465 ± .042) x 10~ 4 degrees/sid. day 2 
(e 0 ) s2 = 7.19 ± .37 km 
( e i)s 2 = (-HI ± .010) km/solar day 
= (.111 ± .010) km/sid. day . 

The mean value of the inclination during the second simulated trajectory period is 


From (70), 


(i ) s2 = 32.836 ± .003 degrees . 


(t 0 ) a = -(14.07 ± .30) days from 10.0 December 1963 . 


From (71), 


(— 0 ) s2 


9.095 


+ .009 
- .009 


degrees west of 50.0 degrees west longitude . 
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From (72), 


(y 0 ) s2 = -(1.300 ± .008) x 10' 3 degrees/solar day 2 

= -(2.257 ± .015) x 10~ 6 rad./sid. day 2 . 

From (74) and the above value of (t 0 ) s2 , 

(a s ) s2 = 42165.6 ± .5 km . 

A graph of these trajectory simulations is seen in Figure A-l. 

Combining the above results for the two simulated trajectories: from (61), 


(A 2 2 ) s 
(A 2 2)^ 


cos 2 (33.005 ± .003) + 1 


cos 2 (32.836 ± .003) + 1 


= [(42165.6 ± .5/42166.3 ± .4) 2 ] 

= .99840 • .00007 

K) s2 - M.x = ["(59.095 ± .009)] - [-(54.782 ± .008)] , 

••• 2 va_ = -(8.626 ± .034) degrees geographic longitude 


Therefore, from (59), the location of the minor equatorial axis with respect to the "synchronous 
longitude" during the first simulated trajectory (54.782 ± .008 degrees west of Greenwich) is 


v\ 


W, 


-tan" 1 / 


sin 

-(8.626 ± .034)] 


(1.300 ± .008) ( 99840 ± -0 0007) - cos 
(1.261 ± .008) 

[-(8.626 ± .034)] 


= 52.5 ± 2.5 degrees east of the minor equatorial axis . 

From (61A), the best estimate of the geographic location of the nearest extension of the equatorial 
minor axis from the simulated trajectory data is 


( y 22) s 


-54.8 - (52.5 ± 2.5) = -(107.3 ± 2.5) degrees geographic longitude . 


This compares well with the input value of (y 22 ) s = -108.0° used to compute the simulated trajec- 
tories. From (60), the best estimate of the triaxial gravity coefficient J 22 from the simulated 
data (according to the theory of this report) is 


(122)5 


| -(2.188 • .013) x 10" 6 


1 

\12tt 2 [sin 2(52.5 ± 2.5)] (6378.388/42166.3 ± A) 2 

cos 2 (33.005 ± .003) + 1" 


2 

J 


= -(1.64 • .03) x 10- 
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The mean equatorial radius used in the simulation is R 0 - = 6378.388 km, the same used to compute 
the "actual" Syncom II orbits from the radar and Minitrack observations. 

The above value of (j 22 ) s compares reasonably well with the input value of (j 22 ) s = -(1.68) 

X 10' 6 used to compute the simulated trajectories. 

The model error implicit in the difference between the reduced and inputed geodetic coeffi- 
cients for the simulated trajectories warrants an adjustment of the J 22 and k 22 reported in Sec- 
tion 7 from the reduction of the "actual" Syncom n orbits. The values below appear sufficient to 
cover all the known uncertainties of this reduction for a triaxial earth: 

j 22 (actual-adjusted) = -(1.70 ± .05) x 10~ 6 

\ 22 (actual-adjusted) = -(19 ± 6) degrees geographic longitude . 

As Appendix B will show, the principal unknown uncertainty of the reduction is the possible in- 
fluence of higher order earth tesseral anomalies on the drift of Syncom n. When all relevant 
higher order anomalies in the earth's gravity potential are evaluated, the adjusted values above 
will probably remain representative for an "average” triaxial potential field sufficient to consider 
for the future design of synchronous satellites. As a guess, the author would increase the upper 
limit of J 22 to about -(1.80) x 10' 6 (based on some of the recent gravity potentials in Appendix B) 
for design purposes, based on an "average" triaxial geoid. A lower limit of j 22 = 1.60 x 10' 6 
for this purpose appears justifiable. The variance in the location of the major equatorial axis for 
the "average" triaxial geoid is not likely to change appreciably from the value quoted for the 
adjusted figure. The author is presently studying these higher order earth- gravity effects. The 
accumulated influence on synchronous satellites of all higher order earth anomalies, is believed 
to be small compared to the 2nd order anomaly. But it appears that close and continuing observa- 
tions on the drift of these satellites will be rewarded in time by revelation of many v of these 
"tesseral" anomalies to about 4th order with an absolute precision almost as good as that reported 
here for the 2nd order effect. 
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Appendix B 


THE EARTH GRAVITY POTENTIAL AND FORCE FIELD USED IN THIS REPORT: 
COMPARISON WITH PREVIOUS INVESTIGATIONS 


The gravity potential used as the basis for the data reduction in this study is the exterior 
potential of the earth derived in reference 10 for geocentric spherical coordinates referenced to 
the earth's spin axis and its center of mass. The infinite series of spherical harmonics is 
truncated after J 44 . The nontesseral harmonic constants J 20 , J 30 and J 40 are derived from ref- 
erence 13. 

The earth radius R 0 used in this study is: 

R 0 = 6378.388 km . 

The earth's gaussian gravity constant used is: 

= 3.9862677 x 10 s km 3 /sec 2 . 

Neither of these values, taken from reference 14, nor the "zonal geoid" of Reference 13, is felt to 
be the most accurate known to date. They are the values used by the GSFC Tracking and Data 
Systems Directorate to calculate the orbit elements of Syncom H from radar and Minitrack ob- 
servations. They were chosen to insure consistency between the data of this study and these pub- 
lished orbits, inasmuch as the "triaxial" reduction for which this study has been undertaken is not 
significantly affected by the probable errors in these values. The second-order tesseral harmonic 
constants used in the simulation studies were 

J 22 = -1.68 x 1(T 6 
\ 22 = -18°. 

These are the values shown on the "tesseral geoid" below (for the j 22 harmonic). At a later point 
in the analysis, the slightly different values reported in the abstract were estimated. The most 
accurate "zonal geoid" is probably that of Kozai (1962) [See Reference 6], with the following earth 
constants; 


R 0 = 6.378.2 km 
^ = 3.98603 x 10 s km 3 /sec 2 . 

The earth's gravity potential (to fourth order, probably sufficient to account for all significant per- 
turbations on a 24-hour satellite) may be illustrated as follows (following Reference 6, Appendix B): 
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The earth- gravity field (per unit test mass) whose potential is (B-l) is given as the gradient of 
(B-l), or 


F = fF r + ^F x + = W E 


+ x *v E 

r d r r cos dk 


1^1 

+ r d<t> 


(B-2) 


or 

F 


r - 'I ' 1 + ( R o/ r ) 2 [ 3 / 2 J 2 o (3sin 2 0-l) + 9J 22 cos 2 <£cos 2 (\-\ 22 ) 

+ 2 (R 0 /r ) J 30 ( 5 sin 2 0-3) (sin<£) + 6(R 0 /r)j 31 (5sin 2 <£~l) cos <p cos (^ _ ^ 31 ) 
+ 60(R o /r) J 32 cos 2 sin <p cos 2(A. - X- 32 j + ®0(®o/ r ) J 33 cos 3 4> cos 3(\ - A. 33 ) 

+ 5/8(R 0 /r) 2 J 40 (35 sin 4 <fi- 30 sin 2 </> + 3) 

+ 25/2(R 0 /r) 2 J 41 (7 sin 2 cos s in c£ cos 

+ 75/2 (R 0 /r) 2 J 42 (7 sin 2 - l) cos 2 0cos2(\-\ 42 ) 

+ 525(R 0 /r) 2 J 43 cos 3 <j> sin 4> cos 3(k-k„) + 525(R 0 /r) 2 J 44 cos 4 <p cos ♦(*-*«)]} 


(B-3) 
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- ^ 2 (R o yr) 2 |f>J 22 cos 4> sin 2 ^ 22 ) + 3/2(R 0 /rj J 31 [5 sin <fi l] sin (/• ~ X. 3 1 ) 

+ 30 (R 0 / r) J 32 cos cf> sincp sin 2 (\ - X 32 ) + 45 (R 0 / r) J 33 cos 2 cp sin 3 (k - \ 33 ) 

+ 5/2 (R 0 /r) 2 J 41 [7 sin 2 <£-3] sin^sin (/-/„,) + 15(R 0 /r) 2 J 42 (7sin 2 ^-l) cos 0 sin 2 (/ - X 4 2 ) 

(B-4) 

+ 315 (R 0 /r) 2 J 43 cos 2 ^sin^sin 3 (\ - ^ 43 ) 

+ 420 (R 0 /r) 2 J 44 cos 3 <p sin 4 - ^ 44 )| * 

= -j(R 0 /r) 2 |-3J 20 sin<^cos<^ + 6J 22 cos^>sin0cos2 (\->v 22 ) 

- 3/2 (R 0 /r) J 30 (5sin 2 </>- l) cos 4> + 3/2 (R„/r) J 31 (l5sin 2 ^- ll) sin cos (k - k 31 ) 

+ 15 (R 0 /r) J 32 (3 sin 2 <p~ l) cos 4> cos 2 (\ - ^ 32 ) 

+ 45 (R 0 /r)j 33 cos 2 ^sin^cos 3 (\ - \ 33 ) - 5/2 (R 0 / r) 2 J 40 (7 s in 2 4> ~ 3) sin^cos 0 

(B-5) 

+ 5/2 (R 0 /r) 2 J 41 (28 sin 4 27 sin 2 0 + 3) cos (*.-*. 41 ) 

+ 30 (R 0 /r) 2 J 42 (7 sin 2 0- 4) cos ^ sin cos 2 (\ - ^ 42 ) 

+ 105 (R 0 /r) 2 J 43 (4 sin 2 </> - l) cos 2 </)Cos 3 (\ - A. 43 ) 

+ 420 (R 0 / r ) 2 J 44 cos3 sint/icos 4 (\ - \ 44 ) j • 


The actual sea-level surface of the earth is to be conceptualized through (B-l) as a sphere of 
radius 6378 km, around which are superimposed the sum of the separate spherical harmonic 
deviations illustrated. To these static gravity deviations, of course, must be added a cen- 
trifugal earth- rotation potential at the earth's surface, to get the true sea level surface (see Ref- 
erence 10). 

From Table B-l and equation (B-3), the fourth-order tesseral geoids of Kaula (September 
1963), Kozai (1962), Izsak (July 1963), and Zhongolovitch (1957) have been evaluated for the longi- 
tudinal perturbation force on a 24-hour satellite with zero inclination at A. = -54.75° over the 
earth's surface (see Table B-2). The harmonics contributing to this perturbation are J 22 , J 31 , 
J 33 , J 42 and J 44 . The results of this comparison are: 
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Tesseral Coefficients in the Earth's Gravity Potential, To Fourth Order, Ab Reported 1915-1963. 
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Table B-2 


Comparison of Longitudinal Perturbation Forces on 
a 24-Hour Satellite, From Five Tesseral Geoids. 



Full Field 
Longitude 
Acceleration 
Units of 

10 -7 g s *J 

Ratio of 
Triaxial (j 22 ) 
Longitude 
Acceleration 
to Full Field 
Longitude 
Acceleration 
(X = -54.75°) 

Zhongolovitch 

(1957) 

7.71 

1.06 

Kozai (1962) 

1.08 

1.28 

Izsak (July 
1963) 

1.27 

1.19 

Kaula (Sept. 
1963) 

1.77 

1.11 

Wagner (this 
reduction: 
March 1964) 

2.21 

? 


*g s is the radial acceleration of earth gravity at the "synchronous” alti- 
tude (g s = 0.735 ft/sec^). 


Judging from the consistency of the 
acceleration ratios among these investi- 
gators, the "actual" value of J 22 (sepa- 
rated from higher order gravity effects) 
is probably somewhat higher than the 
-1.7 x 10~ 6 reported herein. All the geode- 
sists reporting in Table B-2 agree that 
the next most influential earth tesseral 
at "synchronous" altitudes over most of 
the equator is J 33 . 
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Appendix C 


EXPRESSIONS FOR THE INCLINATION FACTOR 


Equation (25) gives the inclination factor in the drift causing tangential perturbation (due to 
equatorial ellipticity) on a 24-hour satellite with a near- circular orbit, as 


F(i) 


AX (max) sin 2 ( i ) 
cos ( i ) + 2 " 


(C-1) 


AX (max.) is the absolute value of the maximum longitude excursion of the figure- 8 ground track of 
the 24-hour satellite (with a near-circular orbit) from the geographic longitude of the nodes. 

From (18), this longitude excursion function is 

AX = X - X Q ~ tan -1 [cos(i)tan$] - 6 (C-2) 

Differentiating (C-2) with respect to the argument angle e, the minimax excursion arguments 
are found from 


d(AX) 

d& 


0 


cos ( i ) sec 2 9 
1 + cos 2 ( i ) tan 2 9 


Solving (C-3) for sin^ at AX(minimax), 


(C— 3) 


sin0 Ai(«n,.«) " [cos(i) + l] 1/2 , from which 


AX(minimax ) 


(C-4) 


(C-4) in (C-2) gives 


AX(minimax) - tan -1 [cos ( i ) sec ( i )] - tan -1 sec (i) . 

Thus, since only the absolute value of ax (minimax) is required, 

AX(max) = tan' 1 [sec(i)j - 45 ° , (C-5) 
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where the tan ' 1 is to be taken in the first quadrant. For example: for i = 30°, (C-5) evaluates 
the maximum excursion as 


A\(max) = 49.1° - 45° = 4.1° . 

The nodal argument angle at this maximum longitude excursion is 

9 [at A/V(max)] = ±49.1° from the nodes! 

The assumption in (20) that the excursion in longitude from the ascending node could be approxi- 
mated by 


AX. - AX(max) sin 26 , 


predicts the maximum excursion argument as 

G [at A\(max)] = ±45° from the nodes. 

This discrepancy in the assumed longitude excursion function is not serious until i >45°, as sim- 
ulated trajectories with variable inclination have borne out. 

(C-5) can be written as 


AX(max) + 45° 
tan [AX(max) + 45° ] 


tan' 1 [sec (i)] , from which 


_. 1 + AX(max) 

sec (i) 1 - AX(max) 


(C-6) 


for i < 45°. Solving (C-6) for Ax (max), 


AX(max) 


1 ~ cos ( i ) 
1 + cos ( i ) 


(C-7) 


approximately. 

Thus the inclination factor becomes approximately 


F( i ) - cos ( 


. . . sin 2 (i) [l - cos ( i )] 

1 ) '2 1+ cos ( i ) 


l-cos(i)] 2 2 cos ( i ) + 1 - 2 cos ( i ) + cos 2 ( i ) 

cos ( i ) + j — - 2 


cos 2 ( i ) + 1 
2 


(C-8) 
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For example: For i = 30°, 


F (i) from(C -i) = .86603 + 4.1/8 x 57.3 = .8750 
F (*)from (C-8) = -8750 . 

The agreement of F(i) from forms (C-l) or (C-8) is good to the third decimal place as long as the 
inclination is less than 45 degrees. At inclinations higher than 45 degrees, however, the drift 
theory following (20) begins to break down because Av(max) is no longer a small angle. 
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Appendix D 


DERIVATION OF THE EXACT ELLIPTIC INTEGRAL OF DRIFT MOTION FOR 
A 24- HOUR SATELLITE WITH A NEAR- CIRCULAR ORBIT: COMPARISON 
OF THE EXACT SOLUTION WITH THE APPROXIMATE SOLUTIONS 
FOR PERIODS VERY CLOSE TO SYNCHRONOUS 


The differential equation (33) of 24 -hour satellite drift is 
analogous to the equation describing the large-angle oscillations 
of a mathematical pendulum (see Reference 11), as in Figure D-l. 

The equation of angular motion of the mass m under the 
constant gravity force mg 0 is 

F t = mg 0 sin# = m(z#) = ml& 

(D-l) can be rewritten as 

8 + (g o /0 sin 6 = 0 . 

From the theory developed in Reference 11 (pp. 327-335), (D-2) 
has an integral 


(D-l) 


(D-2) 



mg 0 


Figure D-l — Configuration of a 
"mathematical pendulum." 


t (time from 9 = 0)= ( i / g 0 ) 1/2 F(k, £) , 


(D-3) 


where F(k, £) is the elliptic integral of the first kind with argument (or amplitude) 
sin 6/2 1 


56 _ sin 


- 1 


sin 6 (max)/2 


j , and modulus k = sin 0 (max)/2. 


Equation (33): 


y + A 22 sin 2y - 0 , 

with maximum libration angle y 0 , can be put in the form of (D-2) by the transformation of the de- 
pendent variable 


6 = 2y , 


(D-4) 
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with the parameter identification 


60/^ 2^22 • 


(D-4A) 


(D-4) implies the identification 

k = sin 7 0 , ^ = sin -1 [sin y/sin y 0 ] . (D-4B) 

The pendulum solution (D-3), under the transformation (D-4) and identifications (D-4A) and (D-4B) 
becomes 

t (time of drift libration from y = 0) 

= (l/2k 22 ) V2 F^siny 0 , [sin -1 sin y/sin y 0 ] j , (D-5) 

♦ 

F(k, ^), in its full integral form, is 

f d<£ 

F ~ Jo (l-k 2 sin 2 si) 1/2 ’ (D-6) 

(where k 2 = sin 2 y 0 , sin 2 <£ = sin 2 y/sin 2 y 0 ) for the drift libration. In particular, when ± = tt/2; 
then y = y 0 ; y 0 = 0, and the pendulum- drift libration has completed a quarter-period. 

Thus, from (D-5) and (D-6), the full period of the long-term drift libration of the 24-hour 
satellite ground track about the nearest minor equatorial axis longitude is 

i -rr/2 

(l - sin 2 y Q sin 2 i/<) 1/2 ' ( D_7 ) 

0 

The adequacy of the Taylor series expansion approximation of the drift motion in #ie neighborhood 
of y 0 , given in equation (51), may be tested against the exact drift solution implicit in (D-5). Table 
D-l below gives the evaluation of F for arguments within 5 0 of y 0 = 60°, using the integral tables 
in Reference 12. 


Table D-l 


Exact and Approximate Drifts of a 24-Hour Satellite from a Stationary Configuration 
60° East of the Earth's Minor Equatorial Axis. 


0 

II 

Gi 

O 

0 

7 

(degrees) 

7' 

(degrees) 

7" 

(degrees) 

4> 

(degrees) 

F 

(rad.) 

AF 

(rad.) 

At 

(days from) 
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60.0 

60.0 

60 

90 

2.1565 

- 

- 


59.0 

59.003 

59.000 

81.7967 

1.8730 

.2835 

41.619 


58.0 

58.014 

58.001 

78.3056 

1.7564 

.4001 

58.737 


57.0 

57.029 

56.999 

75.5595 

1.6671 

.4894 

71.846 


56.0 

56.051 

56.999 

73.1938 

1.5923 

.5642 

82.827 


55.0 

55.077 

54.996 

71.0617 

1.5265 

.6300 

92.487 
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In Table D-l, AF is the change in the elliptic integral from the "stationary” configuration at y = 60° 
or ^ = 90°. At = (l/2A 22 ) 1/2 AF. A 22 was computed from (30A) with the following gravity-earth con- 
stants and for the inclination of Syncom II: 

R 0 = 6378.2 km 
a s = 42166 km 
J 22 = -1.7 x 1<T 6 
i =33°. 

y' gives the drift position as calculated from the first righthand term of (51) alone (the (At) 2 term). 
y" gives the drift position as calculated from the first two righthand terms of (51). The "actual" 
Syncom n drift in mid- August 1963 began, apparently, at a y 0 between 48° and 58° east of the minor 
axis. Thus, the 16 orbits chosen for the first drift period all should be well represented by the 
(At) 2 -only theory, within the rms error of the longitude observations. Similar exact calculations 
for y 0 = 45°, 50°, and 55° confirm the adequacy of the (At) 2 -only theory to apply to the second 
drift-period orbits. They also prove the contention in section 5 that, for reasonably small excur- 
sions from "synchronism," the convergence of the Taylor series (51) is adequate if additional 
terms are included only when they become of a certain minimum significance to the total drift. 
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Appendix E 


BASIC ORBIT DATA USED IN THIS REPORT 


The Orbit elements for Syncom II in Table E- 1 were calculated at the Goddard Space Flight 
Center from radar and Minitrack observations made during the slow drift periods from mid- August 
1963 to February 1964. 

As an example of the estimation of the ascending equator crossing nearest to epoch, consider 
the orbit geometry at epoch (Figure E-l). 

6 January 1964 at 17.0 hours Universal Time (orbit 2-5) . 


Table E-l 


Syncom II Orbital Elements, August 1963 to February 1964. 


Orbit 

# 

Epoch (Universal Time) 
Year-Month -Day -Hour-Min 

Semimajor 

Axis 

(km) 

Eccentricity 

Inclination 

(degrees) 

Mean 

Anomaly 

(degrees) 

Argument 

of 

Perigee 

(degrees) 

Right Ascen- 
sion of the 
Ascending 
Node 
(degrees) 

1-1 

63-8-22-6-12-14 

42164.58 

.00018 

33.083 

24.126 

26.285 

317.569 

1-2 

63-8-26-17-0 

42164.52 

.00016 

33.090 

190.841 

26.099 

317.454 

1-3 

63-8-31-0-0 

42166.02 

.00018 

33.062 

296.125 

30.073 

317.475 

1-4 

63-9-5-0-0 

42166.39 

.00012 

33.064 

333.521 

357.756 

317.362 

1-5 

63-9-9-0-0 

42166.35 

.00015 

33.048 

326.207 

9.077 

317.272 

1-6 

63-9-12-2-0 

42166.55 

.00015 

33.079 

3.657 

4.697 

317.224 

1-7 

63-9-17-2-0 

42166.70 

.00018 

33.043 

12.694 

0.581 

317.165 

1-8 

63-9-20-2-0 

42167.42 

.00018 

33.010 

359.970 

16.282 

317.098 

1-9 

63-9-27-2-0 

42167.51 

.00022 

33.046 

38.922 

344.162 

316.996 

1-10 

63-10-1-2-0 

42168.88 

.00024 

33.039 

26.615 

0.433 

316.944 

1-11 

63-10-8-2-0 

42169.14 

.00020 

33.013 

42.889 

350.866 

316.780 

1-12 

63-10-14-2-0 

42169.78 

.00028 

32.982 

36.727 

2.673 

316.813 

1-13 

63-10-22-2-0 

42171.51 

.00026 

32.993 

62.833 

344.246 

316.603 

1-14 

63-10-30-0-0 

42171.09 

.00028 

32.948 

29.865 

354.548 

316.570 

1-15 

63-11-6-0-0 

42172.15 

.00025 

32.952 

36.699 

354.313 

316.328 

1-16 

63-11-12-0-0 

42172.51 

.00031 

32.920 

108.239 

3.425 

316.308 

2-1 

63-11-28-1-0 

42165.89 

.00005 

32.920 

222.170 

203.901 

315.976 

2-2 

63-12-4-0-0 

42167.20 

.00009 

32.892 

39.435 

17.564 

315.919 

2-3 

63-12-10-0-0 

42167.18 

.00010 

32.881 

51.942 

10.958 

315.877 

2-4 

63-12-16-17-0 

42168.17 

.00007 

32.872 

300.000 

24.505 

315.735 

2-5 

64-1-6-17-0 

42168.01 

.00013 

32.867 

332.997 

11.625 

315.544 

2-6 

64-1-9-6-0 

42169.90 

.00015 

32.857 

165.031 

16.992 

315.469 

2-7 

64-1-20-0-0 

42171.43 

.00012 

32.826 

29.098 

28.842 

315.300 

2-8 

64-1-29-20-0 

42171.91 

.00019 

32.859 

37.956 

13.171 

315.212 

2-9 

64-2-5-16-0 

42172.89 

.00019 

32.800 

321.168 

36.275 

315.075 

2-10 

64-2-10-19-0 

42173.31 

.00014 

32.833 

32.517 

14.553 

314.982 

2-11 

64-2-17-0-0 

42174.89 

.00019 

32.762 

347.774 

35.551 

314.883 
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On 6.0 January 1964, the hour angle of the vernal equinox A west of Greenwich (expressed in hours, 
with 24 hours = 360°) was 

6 hrs 58 m in utes 27.484 seconds (from the Nautical Almanac) 

On 7.0 January 1964, the hour angle of A was 

7 hrs 2 min 24.036 sec. 

Interpolating, the hour angle of A on 6 January at 17 hours Universal Time was 

0 hours 1 minute 15.042 seconds, or 
0.313 degrees west of Greenwich . 

In Figure E-l, the orbit angle 27.003° is taken directly as 360° - the mean anomaly, because the 
orbit is nearly circular. The reported period for this orbit was 

X = 1436.21696 minutes . 

P 

The earth's sidereal rotation period is taken to be 

T e«rth = 1436.06817 minutes . 

Thus, if the satellite is assumed to traverse orbit 2-5 at a nearly uniform rate, it will reach the 
celestial equator at a time when the Greenwich meridian has proceeded eastward from the epoch 

15.378 x 1436.21696/1436.06817 = 15.380° . 

Thus, the estimated geographic longitude of the ascending equator crossing nearest to the epoch 
of orbit 2-5 is 


Ascending equator crossing longitude = -(44.456 + 0.313 + 15.380) 


= -60.149°. 


The estimated time of this crossing is 


15.380°/15° /hr - = 1.025 hours after the epoch . 

The crossing time (Table 1) is thus estimated to be at 

6.751 January 1964 (18.025/24 + 6.0^January 1964) . 
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Appendix F 


CALCULATION OF THE RADIATION PRESSURE ON SYNCOM H 


Consider a flat- plate element AA, of the surface of Syncom H, whose normal n is at angle e 
with respect to the sun's rays (Figure F-l). The sun's radiant energy can be thought of as being 
made up of a stream of material particles such as dm , moving at the speed of light c . If the energy 
of each particle is dE, then from Einstein's classical energy- mass equivalence statement, 

dm = dE/c 2 . (F-l) 

Upon striking the surface, the incident radiation may: 

• Reflect completely off the surface at an equal "reflection" angle, undegraded in energy 

• Be absorbed into the body of the plate as thermal energy and partially reradiated in all di- 
rections from the surface at a reduced flux, depending on the surface and on the thermal 
properties of the plate and body of the spacecraft 

•Be partially "reflected" and absorbed and reradiated, depending on the surface properties of 
the plate 

An estimate of the radiation pressure on Syncom n will be calculated assuming complete light ab- 
sorption with no reradiation. This is not the most conservative condition but will serve to show 
the order of magnitude of the effect. 


SUN 





Figure F-l —Flat plate AA of satel I ite surface, with normal 7 ; 
at angle 0 with respect to sun's rays. 
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Light- particle dm has c dm momentum in the direction of the sun's rays before striking the 
plate element AA. Thus (c) dm momentum is transferred to the plate with each light- particle 
collision. From Newton's second law, the impulse transferred to the plate in the time of action 
dt, of dm alone, is; 


/• dt 

(dF)^ dt = ( c ) d" 1 • (F-2) 

Jo 

dF acts on the plate element in the direction of the sun’s rays. Assume that the discontinuous col- 
lision processes of (F-2) are so frequent as to amount to a smooth transfer of momentum between 
the stream of light particles and the plate; dt can then be replaced by At , a small but finite time 
interval; (dF) din by AF, a smooth, constant small reactive force on the plate element AA; and dm can 
be replaced by Am, a small but finite light particle mass impinging on the plate surface AA in At 
time. (F-2) then becomes 


AF (radiation force) At = (c)Am, or 

Am 

AF(radiation force) = (c) ■ (F-3) 

By the mass- energy equivalence relation (F-l) written for the finite small elements involved in 
the continuous momentum transfer, 


Am - AE/c 2 , (F-4) 

where AE is the energy flux falling on plate element AA in At time. Clearly, 

AE = p AAcos6>(At) , (F-5) 

where p is the sun's energy flux in units of energy /time-area, and AA cos 6 is the projected area 
of the element in the direction of the sun's rays. Combining (F-3), (F-4) and (F-5); 

af/aa = projected area- radiation pressure = p s = ^ ’ (F-6) 

1 in the direction of the sun's rays 

(AA p = AA cos d = projected area of the plate element in the direction of the sun's rays.) 


The value of p outside the earth's atmosphere is estimated to be (see Reference 9) 

P = 95.5 ft-pounds/ft 2 - sec. 

C = 9.835 x 10 8 ft. /sec. 
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(F-6) thus becomes; 


p s (projected area- solar pressure 
against a fully absorbing 
surface in the direction of the 
sun's rays) 


.9725 x 10" 7 pounds/ft. 2 


(F-7) 


Figure F-2 shows the configuration of Syncom II with respect to the sun during the drift. For 
the cylindrical configuration in Figure F-2, from (F-7); 

F (absorb- total) = F (absorb- body) + F(absorb-end) = .9725 x 10' 7 (HD cos 8 + 7tD 2 sin 8 / 4)(F-8) 

The weight of Syncom n in the 24-hour orbit (including the apogee motor) is about 75 pounds. 
Other parameters are: 

H = 15" 

D = 28” 

8 = 21° (in late August 1963) . 



Figure F-2— Configuration of Syncom II with respect to the sun during drift. 
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Thus, HDcos 0 +ttD 2 sin0/4 = 613 in 2 = 4.25 ft 2 ; from (F-8), 


F (radiation force with total absorption) = 4.13 x 10' 7 pounds. 

The mean solar gravity force on Syncom n in orbit is given as 

fi s X 75/32. 15 

F (solar gravity) - * (1 a.u.) 2 ' 

where 1 astronomical unit (A.u. = the earth’s mean distance from the sun) ± 92,900,000 miles 
= 4.9 x 10 n ft. 

M s = 333,0O0 m e = 333,000 x 1.40765 X 10 16 ft 3 /sec 2 = 4.18 x 10 21 ft 3 /sec 2 . 


Therefore 


F (solar gravity) i 4.18 x 10 21 x 75 
on Syncom II (4.9) 2 x 10 22 x 32.15 = 4 ' 07 X 10 2 P ounds • 

It is evident that solar radiation pressure on Syncom II causes perturbations which are insignifi- 
cant compared to solar gravity perturbations. 


56 


